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Abstract 

The triple gauge-boson couplings involving the W are determined using data 
samples collected with the ALEPH detector at mean centre-of-mass energies of 183 GeV 
and 189 GeV, corresponding to integrated luminosities of 57 pb _1 and 174 pb _1 , 
respectively. The couplings, g^ , k 7 and A 7 , are measured using W-pair events, single- 
W production and single-7 production. Each coupling is measured individually with 
the other two couplings fixed at their Standard Model value. Including ALEPH results 
from lower energies, the 95% confidence level intervals for the deviation to the Standard 
Model are 

-0.087 < Agf < 0.141 
-0.200 < Ak 7 < 0.258 
-0.062 < A 7 < 0.147. 

Fits are also presented where two or all three couplings are allowed to vary. In 
addition, W-pair events are used to set limits on the C- or P-violating couplings gj, 
g^, Ry, and Ay, where V denotes either 7 or Z. No deviations from the Standard Model 
expectations are observed. 
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1 Introduction 



The existence of the triple gauge-boson couplings (TGC) is a direct consequence of the 
SU(2) L x U(1) Y gauge theory. The measurement of the TGCs represents a fundamental test 
of the non-Abelian nature of the Standard Model. The triple WW7 and WWZ couplings have 
been studied at LEP in e + e~ collisions at energies above the W-pair production threshold, 
using direct W-pair production (e + e~ — ► W + W~) single-W production (e + e~ — > Weu) 
and single-7 production (e + e~ — * ^7(7)) Measurements of the TGCs have also been 

made at the Tevatron from studies of di-boson production ||. This paper presents new 
results for the TGCs from analyses of W-pair, single-W, and single-7 final states using data 
recorded in 1997 and 1998 with the ALEPH detector. In 1997 and 1998 ALEPH recorded 
total integrated luminosities of 56.81 pb _1 and 174.20 pb _1 , at mean centre-of-mass energies 
of 182.66 GeV and 188.63 GeV, denoted as 183 and 189 GeV. 

The most general Lorentz invariant parametrisation of the WW7 and WWZ vertices 
can be described by 14 independent complex couplings pHl? 7 for each vertex: " x - x 



1 ' 64 ' 

g§ , Ky, Ay, k v and Ay, where V denotes either 7 or Z. Assuming electromagnetic gauge 
invariance, C- and P-conservation, the set of 14 couplings can be reduced to 5 parameters: 
gf, k 7 , Kz, A 7 and Az, with Standard Model values gf = ki = « 7 = 1 and Az = A 7 = 0. 
Precision measurements at the Z resonance at LEP and SLC also provide bounds on the 
couplings [pT0|, [TTj| . However, local SU(2) L x U(l)y gauge invariance reduces the relevance of 



these bounds |TU[ and introduces the constraints: 



Ak z = -Ak 7 tan 2 6 W + Agf , 
Az = A 7 , 

where A denotes the deviation of the respective quantity from its non-zero Standard Model 
value, and #w is the weak mixing angle. Hence, only three parameters remain: Agf, Ak 7 , 
and A 7 ||. 

Using data from e + e~ — > W + W~ final states all three couplings Agf, Ak 7 and A 7 can 
be tested, whereas the single-W and single-7 final states allow measurements of only the 
WW7-couplings, Ak 7 and A 7 . Although the contribution from W-pair production dominates 
the combined limits, the single-W and single-7 events provide complementary information, 
which enhances the sensitivity especially for Ak 7 . 

In this analysis the three couplings Agf, Ak 7 and A 7 are measured individually with the 
two other couplings fixed at zero, their Standard Model value. Fits are also presented, where 
two or all three couplings are allowed to vary. 

The C- or P- violating sector of the TGCs is weakly bound. Indirect limits on k 7 , 
A 7 , kz and gf exist, while there are no direct or indirect limits on the parameters Az, 
gj, gj and gf [p]] . Only the parameter A 7 is tightly constrained by precision low-energy 



measurements |T3|. This paper includes, for the first time, single-parameter fits to the 



unconstrained real and imaginary parts of the 8 couplings g^, g§ , Ky, and Ay, all zero in 
the Standard Model, based on an analysis of semileptonic [evqq and \ivqq) W-pair events. 

The paper is organised as follows. In Section a brief description of the ALEPH detector 
is given. The Monte Carlo event generators used in the analyses are presented in Section ^[ 
The analysis of the single-7 final states is described in Section |] and Section |5| is devoted 
to the single-W analysis. The description of the two analyses is rather concise, as they are 
presented in earlier publications P,ffl. In Section Blthe measurement of TGCs from W-pair 
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events is discussed in detail. Finally, all measurements are combined with ALEPH results 
from W + W~ production at 172 GeV [1[], single-W production at 183 GeV and single-7 
production at 183 GeV ||. The resulting single- and multi-parameter fits are discussed in 
Section [7|, followed by a summary and conclusions in Section |[ 



2 The ALEPH Detector 



A detailed description of the ALEPH detector can be found in Ref. fl3|,|T4||. The central 
part of the ALEPH detector is dedicated to the reconstruction of the trajectories of charged 
particles. Following a charged particle from the interaction point outwards, the trajectory 
is measured by a two- layer silicon strip vertex detector (VDET), a cylindrical drift chamber 
and a large time projection chamber (TPC). The three tracking detectors are immersed in 
a 1.5 T axial field provided by a superconducting solenoidal coil. Altogether they measure 
charged particle momenta with a resolution of 5pr/PT = 6 x 10~ 4 pt © 0.005 (pt in GeV/c). 
In the following, charged particle tracks reconstructed with at least one hit in the VDET, 
at least four hits in the TPC, and originating from within a cylinder of length 20 cm and 2 
cm radius centred on the nominal interaction point and parallel with the beam, are referred 
to as good tracks. 

Photons and electrons are identified in the electromagnetic calorimeter (ECAL), situated 
between the TPC and the coil. It is a lead-proportional-wire sampling calorimeter segmented 
in 0.9° x 0.9° towers read out in three sections in depth. It has a total thickness of 22 radiation 
lengths and yields a relative energy resolution of 0.18/\/E'-|-0.009, with E in GeV, for isolated 
photons. At low polar angles, the ECAL is supplemented by two calorimeters, LCAL and 
SiCAL, principally used to measure the integrated luminosity collected by the experiment. 
Electrons are identified by their transverse and longitudinal shower profiles in ECAL and 
their specific ionisation in the TPC. A detailed description of the photon identification can 
be found in [f[4|| . 

The iron return yoke is equipped with 23 layers of streamer tubes and forms the hadron 
calorimeter (HCAL). The latter provides a relative energy resolution of charged and neutral 
hadrons of 0.85/ y/E, with E in GeV. Muons are distinguished from hadrons by their distinct 
pattern in HCAL and by the muon chambers composed of two double-layers of streamer tubes 
outside HCAL. 

The information from the tracking detectors and the calorimeters are combined in an 
energy flow algorithm fl4]]. For each event, the algorithm provides a set of charged and 



neutral reconstructed particles, called energy flow objects, which are used in the analysis. 
Studies of Z — > qq events show that the angular resolution of jets reconstructed from 
energy flow objects is typically 30 mrad in space and the energy resolution approximately 
o~e = (0-Qy/E + 0.6) (1 + cos 2 8) GeV, where E is the jet energy in GeV and 8 is the polar 
angle with respect to the z-axis along the e~ beam direction. 



3 Monte Carlo generators 

Samples of fully simulated events, reconstructed with the same program as the data, are used 
for the design of the selections, determination of the signal efficiencies and the estimation 
of the background. The size of the generated signal samples correspond to 20 times (for the 
single-7 and single-W analysis) and up to 80 times (for the W-pair analysis) the collected 
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luminosity. 

The efficiency for the single-7 cross section measurement is estimated using a modified 
version of the KORALZ [15| Monte Carlo program. The KORALZ generator simulates initial state 



photons using YFS exponentiation [16[]. The generator is modified to include the effects from 
photons produced as bremsstrahlung off the exchanged virtual W. This treatment includes 
the expected Standard Model contribution and possible anomalous couplings together with 
the interference of the two. The effect on the overall cross section is found to be small 
(~0.2%) for Standard Model couplings. However, it can be as large as a few percent in certain 
kinematical regions. The predictions of the modified KORALZ Monte Carlo are confirmed by 
an independent generator NUNUGPV Jl?]], which is based on exact lowest order amplitudes 
for the production of up to three photons in the final state, modified for higher order QED 
effects using transverse momentum dependent structure functions. 

For the single-W study the GRC4F program |TB[ is used to simulate the four-fermion 
signal process final state evff. The effective QED coupling constant is fixed to be 
«qed = 1/130.2 as suggested in |19|| . For initial state radiation, the photon structure function 
approach is utilised. Final state radiation and tau decays are simulated with PHOTOS p0| 
and TAUOLA [EHI. 



For the analysis of W-pair final states, the KORALW [22] generator, which includes all 
four-fermion diagrams contributing to W + W~-like final states, is used to produce the 
primary reference sample with a W mass of 80.35 GeV/c 2 . The KORALW generator is 



interfaced with JETSET [g3j, PHOTOS gOf, and TAUOLA |]| for fragmentation, final state 
radiation and r decays, respectively. In addition, several samples are generated using the 
double resonant CC03 || diagrams with non-standard values for one coupling at a time, 
to check the reconstruction and TGC determination. Finally, a sample generated with the 
double resonant CC03 diagrams is used to optimise selection efficiencies and parametrise the 
corrections used in the kinematic fitting. 

In order to include the effects from various background processes, Monte Carlo samples 
are generated with a corresponding integrated luminosity of each background sample of at 
least 20 times that of the data. PYTHIA [f23] is used to generate e + e~^ 99(7), ZZ, Zee, 
and e + e~— > Weu event samples. In the ZZ sample, events with W + W~-like final states are 
discarded to avoid double counting. Two-photon processes are simulated with the PH0T02 



generator. The KORALZ [15| and UNIBAB |25[ generators are used for the di-lepton final states. 



4 Single-7 production analysis 

Events with one or more photons and missing energy can be used to probe the anomalous 
WW7 coupling parameters Ak 7 and A 7 . Although the single-7 channel is less sensitive to 
the couplings compared to the W-pair and single-W channels, it provides complementary 
information. A detailed description of the Standard Model processes involved in the reaction 
e + e~ — > 2/177(7) and the modelling of the measured triple gauge-boson couplings can be found 
in 0. 

The sensitivity to the WW7 couplings in the single-7 channel comes from the W-W 
fusion diagram. The Ws exchanged in this t-channel diagram are predominantly at low 
momentum transfer. The single-7 channel is therefore mainly sensitive to Ak 7 because 
contributions from A 7 contain higher powers of the W momenta. Furthermore, the effect 
of anomalous TGCs depends on the energy of the photon. For low energy photons, below 
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the radiative return to the Z peak the sensitivity arises from the interference between the 
Standard Model and the anomalous contribution; this interference has a linear dependence 
on the TGCs. In the region around the radiative return to the Z peak, the sensitivity is 
minimal. For high energy photons, above the radiative return to the Z peak, the dependence 
on the TGCs is quadratic. 



4.1 Event selection and determination of the TGCs 

The events are selected from the 189 GeV data sample using the procedure described in [f26fl . 
In summary, single-7 events are selected by requiring at least one photon candidate with 
# 7 > 20° and p^/E^ eam > 0.1 and no additional activity in form of reconstructed charged 
tracks or energy deposits in the forward regions (below 14°). Events where a photon has 
converted into a e + e~ pair are not considered. 

Anomalous contributions to the WW7 vertex increase the total cross section and lead 
to characteristic energy and angular distributions of the final state photons. For the single- 
7 channel the TGCs are extracted from the data by performing a maximum likelihood fit 
based on the overall number of observed photons, their polar angles # 7 and scaled energies 
x E {= E^/Ebeam) of the form 

log L = log liVexpJ ,. ' + £ log P U (1) 

^obs ! 

where Pi is the probability density function of observing event i with a given value of xe 
and 6* 7 and N eyiV is the expected number of events including background. The probability 
density function and the expected number of events for different values of the couplings are 
constructed by reweighting fully simulated single-7 events. Distributions of the polar angle 
and the scaled energy for single-7 events are shown in Figure |l|. 

Two separate kinematic regions are used in the fit, excluding a region around the Z peak 
return, where the sensitivity is small. Defining Ez — (s — m%)/ (2a/s), the excluded region 
is E^ — 3r z < Ej < £!? + 0.5r z . The total numbers of photons used in the fit are 120 (128 
expected) and 260 (258 expected) below and above the excluded region. 



4.2 Results 

At present energies, the cross section and the shape contribute equally in the likelihood 
function for Ak 7 , whereas the result for A 7 is dominated by the sensitivity to the shape above 
the excluded region. The estimation of the systematic uncertainties follows the procedure 
described in Ref. Q and the different contributions are summarised in Table |1| The fitted 
results for the 189 GeV data for single parameter fits, where each coupling is determined 
setting the other coupling to its Standard Model value, are 

Ak 7 = 0.4 ±0.7 ±0.2 (A 7 = 0) 
A 7 = 0.3 ±0.9 ±0.2 (Ak 7 = 0) 

where the first error is the statistical error and the second is the systematic uncertainty. The 
95% confidence level limits including systematic errors are: 

-1.1 < Ak 7 < 1.8 (A 7 = 0) 
-1.5 < A 7 < 2.0 (Ak 7 = 0). 
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Table 1: Summary of the systematic errors on single parameter fits for Are 7 and A 7 from the single-7 
analysis at 189 GeV. 



bource 


Ak 7 


\ 

A™ 


Acceptance corrections 


0.08 


0.08 


Photon energy calibration 


0.11 


0.14 


Background 


0.05 


0.05 


Luminosity 


0.03 


0.03 


Theoretical uncertainty 


0.13 


0.15 


Total 


0.20 


0.22 



The validity of these 95% CL. limits and the error from the likelihood fit have been checked 
using many Monte Carlo samples corresponding to the data luminosity as described in ||. 

Combining with the previous measurement for centre-of-mass energies between 161 and 
183 GeV [§J, the 95% C.L. limits on Ak 7 and A 7 from single-7 production are 

-1.0 < Ak 7 < 1.5 (A 7 = 0), 
-1.4 < A 7 < 1.8 (Ak 7 = 0). 

The negative log-likelihood functions curves are shown in Figure |2| for the 189 GeV data, 
161 - 183 GeV data M and the combined results. In the combination, the systematic errors 
from acceptance and theoretical prediction are assumed to be fully correlated, while all other 
sources are taken as uncorrelated. 



5 Single-W production analysis 

Single-W production, e + e~ — > Wez/, is sensitive to the WW7 vertex. This sensitivity comes 
from the 7-W fusion diagram, where the momentum transfer is low. As for the single-7 



channel, the single-W channel is therefore mostly sensitive to Ak 7 [27 . 



5.1 Selection 

The analysis of single-W production is performed on the 189 GeV data sample. All W decay 
modes are used and the selection of each W decay final state, described in the following, 
has been optimised for the single-W signal definition used in the previous analysis at lower 
centre-of-mass energies ||: 

6 C < 34 mrad, 

< > 20 GeV and | cos^l < 0.95 for leptonic decays, 
M q q' > 60GeV/c 2 for hadronic decays, 

where 6 C is the polar angle of the scattered electron, and Qi are the energy and polar angle 
of leptons from the W decay. M q q/ is the invariant mass of the quark pair. The cut angle at 
34 mrad corresponds to the lower edge of the acceptance of the ALEPH detector. 

As single-W production is dominated by t-channel processes, the outgoing electron 
is predominantly emitted at small polar angles. Another specific feature is the large 
missing momentum carried away by the electron-neutrino, and therefore a common selection 
criterium for all single-W final states is the requirement of the missing momentum direction 
to be within the detector acceptance, | cos# miss | < 0.9. 
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5.1.1 Leptonic Selection 

The leptonic W decay is characterised by a high energy isolated lepton. Allowing for a multi- 
prong decay of the tau, events with one or three good tracks (| cos6*| < 0.95) are accepted. 
The selection cuts are the same as in the analysis in Ref. || and are summarised here. 

In addition to the cut on the missing momentum direction, tagged two-photon events are 
rejected by requiring that no energy be detected within a cone of 12° around the beam axis 
(Eia = 0). 

The remaining backgrounds, mainly untagged two-photon events and two-fermion events, 
are eliminated by requiring that the transverse missing momentum be greater than 0.06i/i. 
This threshold is increased to O.Ia/s if the missing momentum direction points to within 10° 
in azimuth to the boundaries between the two LCAL halves or between the six inner sectors 
for the TPC. It is required that no energy is found within a wedge of 10° opposite to the 
direction of the lepton transverse momentum. To reduce the background from e + e~ — > Zee 
with Z decaying to neutrinos, events are rejected if an electron candidate track is identified 
and its energy, including the neutral energy in a 10° cone around the track, is less than 
20 GeV. 

The selection efficiencies for the three final states are 75% [evev\ 77% (ez//iz/) and 43% 
(ez/rz/). The main background source is Zee where the Z decays to z/ M z/ M or z/ T z/ r (the z/ e z/ e 
case is a four-fermion final state which is e + e~ — > Weu like and is part of the signal). 

In the data, 23 events are observed in agreement with the expectation from the Standard 
Model of 26.5 events (17.7 signal events). The composition is 15 events with an electron 
(8.4 signal and 6.2 background expected), 4 events with a muon (6.6 signal and 0.4 
background expected) and 4 events with a tau (2.7 signal and 2.2 background expected). 
The distributions of the lepton transverse momentum, p l T , for single- W events passing the 
final selection cuts for the leptonic electron and muon W decay, are shown in Figure |3|. 

5.1.2 Hadronic selection 

For the hadronic W decay, the event topology is characterised by two acoplanar jets with an 
invariant mass around that of the W boson. The selection is the same as in the analysis in 
Ref. 0. In addition to the cut on missing momentum direction, at least seven good tracks 
are required. Similarly to the leptonic selection, tagged two-photon events and two-fermion 
events with initial state radiation are rejected by demanding that the energy E 12 be less than 
0.025-y/s. The visible mass is required to exceed 60GeV/c 2 and to be less than 90 GeV/c 2 
to reject untagged two-photon events at the low end of the mass spectrum and ZZ events at 
the high end. 

Events for which the energy in a wedge of 30° centred on the transverse missing 
momentum direction is greater than Cly^ are rejected. The acollinearity angle between 
the two hemisphere (defined by the event thrust axis) momentum directions is required to 
be less than 165°. 

The semileptonic final state [tvqq) of W-pair production is efficiently rejected by requiring 
that no identified electron or muon with an energy of more than 0.05^/s be reconstructed. 
The tau jet reconstruction algorithm of Ref. || is used in order to further reject semileptonic 
decays of W-pairs which contain a tau lepton. After all cuts, the semileptonic W-pair 
production, primarily events with one tau lepton, remains as the dominant background. 

The efficiency for the hadronic W channel is about 43%. In the data, 53 events are 
observed, in agreement with the Standard Model expectation of 63.1 events (23.5 signal 
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Table 2: Summary of the relative systematic uncertainties in the expected numbers of selected signal events 
of the leptonic and the hadronic channels from the single- W analysis at 189 GeV. 



bource 


a a T exp i A T exp 


a a T exp i AT exp 
AA had /A had 


Luminosity 


±0.01 


±0.03 


Calorimeter calibration 




+0.11 
-0.08 


E 12 inefficiency 


±0.01 




Signal and background cross-section 


±0.06 


±0.05 


Fragmentation 




±0.05 


Total 


±0.07 


+0.13 

-0.11 



events). The visible mass distribution of the selected events is displayed in Figure 0. 
5.2 Results 

Limits on Ak 7 and A 7 are derived from the total rate of single-W events, which is sensitive 
to the WW7 couplings. The upper limit on the single-W signal cross section has been 
calculated while varying only one coupling at a time, and the 95% CL. limits on Ak 7 and 
A 7 for the 189 GeV data are 

-2.09 < Ak 7 < 0.20 (A 7 = 0) 
-0.77 < A 7 < 0.79 (Ak 7 = 0), 

including the systematic uncertainties. The different contributions to the systematic errors 
are summarised in Table |2|. The total systematic error amounts to 7% for the leptonic 
and 12% for the hadronic channel on the predicted numbers of signal events. The overall 
systematic errors are small compared to the statistical precision, which amounts to 33% for 
the leptonic channel and 55% for the hadronic channel. 

Combining with the previous measurement for centre-of-mass energies between 161 and 
183 GeV [f|, the 95% C.L. limits on Ak 7 and A 7 from single-W are 

-2.12 < Ak 7 < 0.23 (A 7 = 0), 
-0.76 < A 7 < 0.78 (Ak 7 = 0). 

The corresponding logL curves are shown in Figure |] for Ak 7 (A 7 = 0) and A 7 (Ak 7 = 0) for 
the 189 GeV data, 161 - 183 GeV data and the combined results. 

6 W-Pair production analysis 

The large number of W-pair events produced yield the dominant sensitivity to the TGCs. 
The process is sensitive to both the WW7 and the WWZ couplings via the s-channel W- 
pair production diagrams and the sensitivity to the coupling A 7 is higher than that of the 
single-W and the single-7 processes. 

6.1 Event selection and kinematic reconstruction 

In this section the event selections for the three distinct W + W~ event topologies, iuqq, 
qqqq, and £i/£v, are described. Selected events are exclusively classified in the following 
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Table 3: The numbers of events after all cuts applied in the final W + W _ TGC results for data and Monte 
Carlo simulation in all channels at centre-of-mass energies of 183 and 189 GeV. The number of Monte Carlo 
events is normalised to the respective integrated luminosity of the data. The quoted efficiencies e and purities 
p are determined from CC03 events with mw = 80.35 GeV/c 2 . For a given W + W _ channel, contributions 
from other channels arc considered as background. 





y/s = 183 GeV 


y/s = 189 GeV 




evqq 


fivqq 


rvqq 


qqqq 


tvtv 


evqq 


\ivqq 


rvqq 


qqqq 


ivtv 


N 

11 exp 
Ndata 


105.6 
117 


107.1 

95 


90.5 

88 


296.8 
314 


29.7 
29 


360.6 
361 


369.9 
370 


230.2 
224 


1202.9 
1130 


101.3 
102 


e(%) 
P(%) 


76.2 
96.3 


79.5 
97.8 


50.7 
73.5 


67.0 
90.0 


61.9 
89.8 


74.0 
96.7 


78.1 
98.1 


44.8 
80.6 


78.6 
86.2 


61.9 
90.3 



order of priority: fMvqq, evqq, QQQQ, TV QQ-, and Ivlv. The expected numbers of events after 
all cuts used in the TGC results for signal and background processes at both centre-of-mass 
energies are summarised in Table [3] for each channel, along with the corresponding selection 
efficiencies and purities. 

6.1.1 W ~ W — > evqq and W + W _ — > n>vqq events 

The event selection procedure for semileptonic W + W~ events is similar to that used for the W 
mass measurement at the corresponding energy P8"[2"9||. At 183 GeV, events are reconstructed 



such that they contain a high energy lepton candidate and two jets . The charged particle 
with the highest momentum component anti-parallel to the missing momentum is chosen as 
lepton candidate. At 189 GeV the selection criteria for the lepton track are slightly changed, 
using the lepton track isolation |2"9]]. The DURHAM-PE [EST] clustering algorithm is applied to 



all energy flow objects not used to construct the lepton four-momentum, and these are forced 
into two jets. After this preselection, the probability for the event being signal is determined 
using the momentum of the lepton, the total missing transverse momentum and the lepton 
isolation from the closest jet. 

At this stage events passing a cut on the probability are subjected to a kinematic fit in 
order to improve the resolution on the reconstructed four-momenta of the W decay products. 
The kinematic fit and additional reconstruction cuts are described in the context of the 
specific TGC analysis. 

A W-pair event can be characterised by the five measured angles, 9w, the W~ production 
angle between the W~ and initial e~ in the W + W~ rest frame, the polar and azimuthal angles 
of the lepton, 9* and </>*, in the rest frame of its parent W and the polar and azimuthal angles 
of a quark jet, 9? et and 0*j et , m the rest frame of its parent W. The distributions of the five 
angles cos 9w, cos 9*, <f>*, cos 9* et and 0* et , for evqq and \ivqq events at 189 GeV after selection 
and reconstruction is represented in Figure |5]. 

6.1.2 W+W - — >ruqq events 

The event selection is based on two complementary algorithms, developed for the cross 
section measurement at 161 |5U| and 172 GeV |3~I] , but modified to account for the change in 



event kinematics with centre-of-mass energy. In summary, events passing a set of preselection 
cuts are selected as semileptonic r candidate events if they fulfil either a global or topological 
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selection. The DURHAM-PE clustering algorithm is then applied to all energy flow objects that 
are not used to construct the tau four-momentum, and these are forced into two jets. More 
detailed descriptions of the selection and tau reconstruction can be found in the publications 



on the W mass determination at the corresponding energies [p8, 29 . 

To improve the resolution of the angular observables a 3-constraint kinematic fit is 
applied, requiring four-momentum conservation and reference mass constraints. In the 
kinematic fit the direction of the r is approximated by its visible decay products and the extra 
energy loss is compensated by correction coefficients obtained from Monte Carlo simulated 
W + W~^rz/gg events. For single prong r decays the charge of the r is directly accessible, 
but in the case of three-prong r decays ambiguities arise due to mis-assigned particles from 
the jets to the r. For three-prong r decays the charge of the r is therefore determined 
from the sign of the pseudorapidity-weighted average jet charge of the r decay products 



(Section |6.1.3| ), where the pseudo-rapidity is defined with respect to jet-axis. The charge 
mis- assignment in rvqq events is 5% for one-prong and 41% for three-prong r decays. 

The distribution of the cosine of the W~ production angle from rvqq events can be seen 
in Figure |. 

6.1.3 W+W - — >qqqq events 

To extract the hadronic W + W~ signal with high purity and efficiency, the selection is 
based on a neural network [|32]. Events passing a preselection designed to remove the qq(7) 
background, are assigned a neural network output, based on global event properties, heavy 
quark flavour tagging, jet properties and WW kinematics. A detailed description of the 



selection algorithm at 183 and 189 GeV is given in p2|. 

For the hadronic W + W~ events the reconstruction of the relevant information is more 
complicated since the W~ direction is not known and the information on the particle flavours 
in either W system is not discriminant. In this case the four jets can be paired in three 
different ways. To select the best pairing, a 6-constraint kinematic fit is applied to all three 
possible pairings. The kinematic fit requires four-momentum conservation and reference 
mass constraints. The four-momenta obtained in the kinematic fit for the pairing with 
the lowest x 2 value are then used in the final determination of the TGCs, while the other 
combinations are discarded. The efficiency of this algorithm to find the correct combination 
was found to be 78% at 183 and 75% at 189 GeV. 

To assign a jet pair to the W + or W~ a jet charge algorithm is used. The jet charge, 
Qjet, is obtained from the pseudorapidity-weighted average charge of jet particles. The jet 
pair charge is defined by the sum of the two jets assigned to a W, Qw = Qjeti + Qjet2- A jet 
pair is then assigned to the W + based on the charge difference between the two jet pairs, 
AQ, with a probability P + . The probability P + is given by 

P (AO) - N + (AQ) 

P + {AQ)- Af+{AQ)+Af+{ _ AQ y (2) 

where A/+ is the probability density function for the charge difference between the two W 
systems for true W + jet pairs obtained from Monte Carlo event samples 0. Figure [7] shows 
the distribution of the di-jet charge of the two W systems for true W + and W~ decays, 
obtained from Monte Carlo generator information. The distribution of the di-jet charge 
from semileptonic events is also shown for both data and Monte Carlo. The data are well 
reproduced by the Monte Carlo simulation. The charge assignment efficiency for correctly 
paired hadronic W + W~ events amounts to approximately 76% for P + > 0.5. 
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The distribution of the cosine of the W production angle from qqqq events can be found 
in Figure []. 



6.1.4 W+W~^£v£v events 

The selection of W + W~— *iviv events (where I denotes an electron or muon) is mainly based 
on variables used in Ref. p0| , namely missing transverse momentum, missing mass and 



kinematic properties of the lepton candidate. The information of these and other variables 
is combined in a neural network. A detailed description of all the variables used in the neural 
network can be found in the appendix Kj. 

For purely leptonic W + W~ events the momenta of the two neutrinos are unknown. 
However, in the absence of ISR and neglecting the W width, the constraint that the two 
Iv systems should have the W mass (Mi = M2 = 80.35 GeV/c 2 ) in combination with 
the usual four-momentum conservation allows a reconstruction of the neutrino momenta. 
The quadratic nature of the mass constraint results in a two-fold ambiguity, corresponding 
to flipping both neutrinos with respect to the plane defined by the charged leptons. As 
detector resolution, ISR and the finite W width are not included in this reconstruction 
hypothesis, 28.0% of the events have no physical solution and a zero-constrained kinematic 
fit is employed. The fit determines a set of values for the reference masses, Mi and M 2 , 
in the mass-constraints for the two Iv systems, for which a physical solution exists. By 
this method 92% of the events without a solution are recovered, resulting in a 97.7% 
reconstruction efficiency for signal events. The majority of the events which fail to have 
a solution are purely leptonic W + W~ events with at least one leptonically decaying r, which 
is the dominant background. In addition to the selection by a cut on the neural network 
output, purely leptonic events are only accepted if a physical solution is found inside the 
mass window of 55 — 105 GeV/c 2 . 

The distribution of the cosine of the W~ production angle for Ivlv events at 189 GeV 
after selection and reconstruction is represented in Figure 19. 



6.2 Determination of the TGCs 

Three different methods, described in the following, are used to extract the couplings in the 
different W-pair final states. 



6.2.1 The optimal observable methods 

The general idea of optimal observables (00) [^] is to project the sensitive kinematic 
information for a given coupling onto the one-dimensional distribution of a suitably defined 
variable O^. The coupling $ can then be extracted from a fit to this distribution or, 
equivalently, from the measurement of the mean value (O^ ) of the optimal observable. 

Since the amplitudes are linear in the TGCs the differential cross section can be expanded 
in these couplings g^ containing no terms beyond the second order 

% = s„(n)(i + E oP -9i + Yl off ■ *9i)> ( 3 ) 

i ij 

where gi denotes any type of couplings and Q denotes phase space variables taking into 
account reconstruction ambiguities for the individual W + W~ channels. The zero-order term, 
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Sq(Q), represents the Standard Model contribution. Using the first order term, a given set 
of couplings, g, can be determined by minimising 

X 2 (9) = E«^ (1) > " E[0V\)V{0)7l{(0f) - E[Of}) } (4) 

u 

where (O^) and V(0) are the measured mean values and their covariance matrix. The 
expected mean values, E[0\ ], are obtained by reweighting of fully simulated W + W~ events. 
In order to ensure maximal sensitivity, gi can be determined by two different approaches: 

• An iterative procedure (denoted OOi in the further text), where the cross section in 
Eq. H is expanded about a given coupling value and consequently the definition of the 
observable Oj is re-optimised. 

• Adding the information contained in the second order term of the expansion in Eq. ^| 
(called OO2 in the further text). This is achieved by including the second optimal 
observable and adding additional terms of the same structure as the ones in Eq. |4] to 
the x 2 j including new terms describing the correlation between and . The 
second order observable increases the sensitivity when the information contained in 
the first order observable decreases p4 |. The covariance matrix for the mean values 



are obtained by reweighting fully simulated W + W events. 

In both cases the information from the measured cross section is included by adding a Poisson 
term to the likelihood function. 

By construction, these methods are bias-free and take into account any experimental 
effect, provided that the Monte Carlo simulation describes the data correctly. For a 
given channel, contributions from any other channel are considered as coupling dependent 
backgrounds. In addition the efficiencies and purities of each selection are parameterised as 
function of the couplings. 

For semileptonic events, both 00 analyses, OOi and OO2, apply a two-constraint 
kinematic fit using four-momentum conservation, equal W-mass hypothesis and including 
a massless neutrino. The corresponding x 2 -probability of the fit is required to be larger 
than 0.005 for an event to be selected for the TGC extraction. This cut improves the purity 
of the sample and discards poorly reconstructed events. In addition, an OO-window cut is 
applied, optimised to improve the sensitivity to the TGCs by reducing the contamination of 
background events with OO values incompatible with W-pair production. 

For hadronically decaying W's, there remains a twofold ambiguity since the quark flavours 
are undetermined. Hence, for semileptonic euqq, [ivqq and Tvqq events the contributions are 
averaged over the quark and anti-quark directions. For hadronic events, considering the 
W charge as undetermined, there is an eightfold ambiguity. This is included in the final 
extraction of the TGCs, where each contribution is weighted with the corresponding di-jet 
charge probability, Eq. [| For leptonic events, iviv, the contributions are averaged over the 
two solutions for the neutrino momenta. 



6.2.2 Maximum likelihood- fit 

A maximum likelihood analysis (LL), Eq. [I], of the channels evqq and fiuqq is used to measure 
the C- or P-violating couplings and as a cross-check for the CP-conserving couplings. The 
measured variables are the five angles described in Sect. |6.1.1| . As in the optimal observable 
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Table 4: The numbers of events after all cuts applied in the maximum likelihood TGC analysis for data 
and Monte Carlo simulation at centre-of-mass energies of 183 and 189 GeV. The number of Monte Carlo 
events is normalised to the respective integrated luminosity of the data. The quoted efficiencies e and purities 
p are determined from CC03 events with mw = 80.35 GeV/c 2 . Only non-W + W~ events are considered as 
background in the calculation of the efficiencies e and purities p. 





y/s = 183 GeV 


y/s = 189 GeV 


Channel 


evqq 


fiuqq 


evqq 


fiuqq 


N 

1,1 exp 
Ndata 


91.8 
98 


97.4 
86 


293.0 
275 


311.7 
310 


e(%) 

P(%) 


66.4 
98.5 


71.0 
99.4 


66.1 
98.6 


70.9 
99.5 



analyses, no quark flavour tagging is performed and the quark and anti-quark directions are 
averaged. The probability density function, P, is given by 

m , _ m±*M , 5) 

where g denotes a set of couplings and the angles fl = (6^y, #f , 0f , #* et , 0* et ) are calculated 
using the charged lepton, neutrino, and quark jet four-vectors. The quantity b(fl) is the 
background distribution as predicted by Monte Carlo. The W + W~ signal distribution, 
s(D,,g), is defined by 

s(ti,g) = J ^^ tru V(0,^ rue )e(^^ true )F(s)^(s,O true ^), (6) 

where s is the squared invariant mass of the W + W~ system, r(f2,f2 true ) is the detector 
resolution function, e(s, f2 true ) is the detection efficiency, F(s) is an initial state radiation 
function |j35| |, and ^(s, f2 true , g) is the lowest-order narrow-width differential cross-section 



for W + W~ production and decay 0. 

The normalisation factors B and S(g) are the integrals of b(Q) and s(Q, g), where S(g) is 
evaluated by reweighting W + W~ Monte Carlo events in order to include detector resolution 
and efficiency. The proper evaluation of the normalisation S(g) is crucial for the success of 
the likelihood method. 

Some approximations to s(Q,g) are made when evaluating the numerator of Eq. [5[ In 
particular, most of the information in the detection efficiency function e(s, 0) arises from 
the charged lepton momentum pi and polar angle d\ in the laboratory reference frame. This 
dependence has been parameterised with a two-dimensional efficiency function. 

The event selection for this method is the same as for the optimal observables. A 
kinematic fit using the equal mass hypothesis and four-momentum conservation is applied, 
where the corresponding x 2 -probability is required to be larger than 0.02. The efficiencies, 
purities and numbers of events for data and Monte Carlo simulation are shown in Table 

6.3 Results 

6.3.1 Comparison of methods 

The three methods (OOi, OO2 and the maximum likelihood) discussed above are used to 
determine the TGCs from WW events. 
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Table 5: Comparison of the three couplings Agf , Ak 7 and A 7 , using the three methods OOi, OO2 and 
maximum likelihood. The error intervals for each coupling are statistical only. 







Method 


Channel 


Coupling 


OOi 


oo 2 


LL 




Agf 


nn7 +o.n 
u - u '-o.io 


10+0- 09 

u - iU -o.io 


O03+O.IO 
u - Ud -0 .10 


evqq 


Ak 7 


—0.34 


0.45+0-35 
—0.33 


0.11+0-47 
—0.30 




A 7 


04 +0.12 
u - u -0.11 


22+0.11 
u - -0.10 


011+041 




Agf 


05+0+0 

u - u ^-o.io 


ooo+o- 10 

u ' uu -0.08 


06+0- 09 
u - UD -0.09 


fivqq 


Ak 7 


-0.02+S-q5 
—0.34 


0.24+0-59 
—0.35 


0-38+0-49 
—0.38 




A 7 


-0 03+ - 10 


-0 08+0- 09 
u - U8 -0.08 


-0 08+0- 09 
u - U8 -0.09 




Agi 


0^1+0.19 
U - Oi -0.29 






rvqq 


Ak 7 
A 7 


-0.71+0-39 
—0.32 

u - uu -0.14 








Agf 


Ofi+0- 10 
U ' Ub -0.09 


_ n 03+0.11 

u - Ud -o.io 


- 


qqqq 


Ak 7 


-0.11+0-30 
—0.27 


0.21+0-34 
—0.78 






A 7 


15+0.11 

u - i5 -o.io 


02+0- 14 
U ' U ^-0.12 






Agi 


17 +0.30 
u - i '-0.20 






tvtv 


Ak 7 
A 7 


_ 35+O.8O 
u - dO -0.41 

0^+0-13 
U - U5 -0.13 







The OOi method is used to measure the couplings Ag 7 , A« 7 and A 7 for all five WW final 
states considered in this analysis, namely evqq, fivqq, Tvqq, qqqq and tvtv. The 00 2 method 
is employed to measure the couplings Ag^ , Ak 7 and A 7 for the evqq, \ivqq and qqqq channel. 
The maximum likelihood method is used in the two semileptonic channels, evqq and [ivqq, 
to measure the standard set of couplings Agf , Ak 7 and A 7 and the real and imaginary parts 
of the C- or P-violating couplings gj , g^, ky, and Ay, where V denotes either 7 or Z. 

A comparison of the three couplings Ag-f, Ak 7 and A 7 , using the three methods OOi, 
OO2 and maximum likelihood, is given in Table [5]. For OOi and OO2, the combined results 
at 183 and 189 GeV are extracted by adding up the corresponding \ 2 terms, while for the 
maximum likelihood method results are extracted by summing up the corresponding logL 
functions. The error intervals for each coupling are defined as the 68% confidence intervals 
obtained by integration of the likelihood functions, to accommodate cases with non-parabolic 
behaviour of the log-likelihood function. 

The linearity of the three fitting procedures is checked by repeating the fits using Monte 



13 



Table 6: The expected error for the three methods, OOi, OO2 and maximum likelihood, at 189 GeV for 
the three couplings Agf, Ak 7 and A 7 . The maximum likelihood method has only been applied in the two 
semileptonic channels, evqq and [ivqq. 







Expected error 


Channel 


Coupling 


OOi 


oo 2 


LL 




Agf 


0.11 


0.12 


0.11 


evqq 


Ak 7 


0.48 


0.43 


0.40 




A 7 


0.13 


0.12 


0.13 




Agf 


0.11 


0.11 


0.10 


\ivqq 


Ak 7 


0.48 


0.41 


0.37 




A 7 


0.11 


0.11 


0.13 



Carlo event samples generated with non-zero values for the TGCs. In all cases the results 
are consistent to those values within the statistical uncertainty of the Monte Carlo samples. 

The consistency of the Monte Carlo simulations with the data is verified by comparisons 
of the distributions of the input quantities to the selections for data and Monte Carlo. In 
addition, the stability of the analysis with respect to the event selection is tested by varying 
the main selection criteria within reasonable limits. In no case significant discrepancies are 
found. 

The reliability of the errors from each fitting procedure is investigated by performing fits 
to a large number of independent Monte Carlo samples, each corresponding to the integrated 
luminosity of the data. These samples, typically 300, are then reweighted to non-zero values 
for the TGCs and passed through the analysis chain. The expected 68% confidence levels, 
obtained from the distributions of the fit values, show good correspondence with the 68% 
confidence intervals obtained for data. 

The consistency of the data results from the three different methods used for the evqq 
and fivqq channels, has been checked by performing fits to a large number of independent 
Monte Carlo samples. The results obtained with the three different methods are compatible 
and the expected spread between the methods show a good agreement with the observed 
differences in the data. 

The expected 68% confidence level intervals, obtained from the distributions of the fit 
values for the three methods at 189 GeV for the three couplings Ag^, Ak 7 and A 7 , are listed 
in Table ^. The expected errors of the OO2 and the maximum likelihood method are very 
similar, whereas the OOi method is slightly worse in the case of Ak 7 . For the final results, 
the OO2 method is therefore employed in the analysis of the evqq, \ivqq and qqqq final states 
and the OOi analysis is used for the remaining rvqq and £v£v final states. 

6.3.2 Systematic uncertainties 

In the following the different sources of systematic errors for each decay channel and their 
determination for the various methods are briefly described. The different contributions 
of each source to the total systematic error for the three couplings Ag^, Ak 7 and A 7 , as 
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obtained with the 00 methods, are given for each channel in Table |7]. The systematic 
uncertainties for the combined evqq and [ivqq channels are listed in Table || for the C- or 
P-violating couplings, determined with the maximum likelihood method. 

The following sources, listed in the approximate relative importance, have been 
considered to be fully correlated between the channels: 

- Fragmentation: The effect of fragmentation in hadronic W decays is estimated by 
fitting the couplings in samples of events generated with KORALW where the default 
JETSET fragmentation is replaced by HERWIG 0] fragmentation. The HERWIG 



fragmentation parameters are tuned at the Z using hadronic events with flavour 



tagging |]29 



W + W cross section: The uncertainty due to the theoretical error on the expected 
W + W~ cross section predicted by KORALW is estimated by changing the W + W~ cross 
section by ±2% [|. 

Luminosity: The effect of the error on the integrated luminosity is estimated by 
varying the measured integrated luminosity by ±0.7% |37[. 



LEP energy: The uncertainty on the LEP energy affects the determination of the 
couplings via the kinematic fitting procedure and the cross section measurement. The 



values of the LEP energy are varied in the range ±0.050 GeV [38|, which has a negligible 
effect on the results. 

W mass: The analysis is repeated using Monte Carlo samples generated at different 
values of the W mass to investigate the effects due to the uncertainty of ±62 MeV in 



the W mass measured at hadron colliders 39 



Calorimeter absolute scale: The absolute energy scale of the electromagnetic and 
hadronic calorimeters is determined using hadronic Z events. The uncertainties in the 
absolute scale are found to be ±0.9% and ±2% for the electromagnetic and hadronic 
calorimeter, respectively. The effect of a possible miscalibration of the calorimeters is 
evaluated on Monte Carlo samples by scaling the electromagnetic and hadronic part 
of the measured energy independently by these amounts. The largest of the observed 
shifts for each calorimeter is combined in quadrature. 

Particle tracking: The definition of a good charged track has been tightened in 
this analysis in order to minimise possible effects from residual tracking distortions 
primarily in the forward regions of the detector. Corrections for the distortions 
are determined by studying Z— > ji~ events, and possible remaining distortions 
have been estimated using Bhabha events. The systematic uncertainty related to 
tracking is assessed by applying the corrections and adding the remaining distortions 
independently to Monte Carlo event samples and repeating the analysis |28[ . 



Jet energy corrections: Detailed comparisons of reconstructed jets in Monte Carlo 
and data are used to parametrise small corrections to Monte Carlo jet energies as 



function of the jet polar angle to the beam axis ||28|| . In order to evaluate the effect from 
the uncertainty in the Monte Carlo jet energy corrections, two alternative correction 
functions, corresponding to ±lcx errors of the discrepancies, are used. The largest shift 
with respect to the nominal correction is taken as the systematic error [EM]. 
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Higher order terms: The effect from missing higher order terms, 0(a 3 ), in the 
simulation of initial state radiation in the KORALW generator, is assessed following the 



procedure described in Ref. |[28|| . In summary, the error on the couplings is determined 
by comparing fits of Monte Carlo samples with events weighted to 0(a)/0(a 2 ) with 
fits to the corresponding unweighted samples, which have been generated in the second 
order leading-log approximation. Recently, new improved calculations with the Double 



Pole Approximation jMj have appeared. The improved CC03 cross section and angular 



distributions predicted by two independent Monte Carlo programs, RacoonWW |yj and 



YFSWW |42|], could introduce small changes on the couplings. For the time being these 



effects have not been included. 
Errors assumed to be uncorrelated between channels include: 

- Bose-Einstein correlations: The effect of Bose-Einstein correlations in the qqqq 
channel is investigated by repeating the analysis on Monte Carlo events generated with 
KORALW and fragmented using JETSET with Bose-Einstein correlations for all particles, 
following the implementation in LUBOEI [0 . The scheme for restoring four-momentum 



conservation denoted BE3, which has been tuned to the LEP1 Z data, is considered [29|. 



Colour reconnection: The uncertainty arising from possible colour reconnection 
effects is assessed by studying Monte Carlo implementations of different colour 



reconnection scenarios in the parton evolution scheme in JETSET [44]. The analysis 
is repeated with qqqq events generated with the EXCALIBUR generator and hadronised 
with and without colour reconnection in the model referred to as SKI, as described 
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2HJ . The systematic error is taken as the difference in fitted couplings from samples 
without colour reconnection and with colour reconnection in about 30% of the events. 

- Background estimation: The error on the couplings from the uncertainties in 
the background estimation is evaluated by varying the normalisation of the main 
background processes. The background from QCD is changed by ±5% based on 
comparisons between data and Monte Carlo simulation. The background from 77, 
Zee and ZZ processes is varied by ±30% ±20% and ±2%, respectively, to account for 
the theoretical uncertainty in the prediction for those processes ||. 

- Monte Carlo statistics: The effect of the Monte Carlo statistics is included in the 
systematic uncertainty. 

- Jet charge assignment: To investigate the effects from the uncertainties on the 
jet charge, the reconstructed W charge is shifted by 0.01. This number is based on 
comparisons between data and Monte Carlo simulation (Z peak data) El. The effect 
on the couplings is found to be negligible. 

The systematic uncertainties listed above as fully correlated between channels are also 
assumed to be fully correlated between years. In addition, the systematic errors from Bose- 
Einstein correlations and colour reconnection are taken to be to fully correlated between 
years. 

For both optimal observable methods, OOi and OO2, the systematic errors have been 
calculated based on the changes in the mean values of the respective observables. This is 
incorporated in the TGC extraction by including the corresponding covariance matrix for 
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Ag? 


Ak 7 




Source 


evqq 


ixvqq 


Tvqq 


qqqq 


ivlv 


evqq 


livqq 


Tvqq 


qqqq 


Ivlv 


evqq 


livqq 


Tvqq 


qqqq 


Ivlv 


Correlated errors 




Fragmentation 


0.01 


0.01 


0.09 


0.02 




0.06 


0.04 


0.20 


0.05 




0.01 




0.15 


0.04 




W + W~ cross section 




0.02 


0.10 


0.02 


0.03 




0.15 


0.09 


0.04 


0.07 




0.02 




0.03 




Luminosity 






0.05 




0.01 




0.03 


0.04 


0.03 


0.03 








0.01 


0.01 


LEP energy 


















0.01 














W mass 






0.03 




0.02 


0.02 




0.17 


0.02 


0.10 






0.06 




0.03 


Calorimeter scale 






0.12 




0.03 


0.03 




0.12 




0.04 


0.01 




0.03 




0.03 


Tracking 










0.04 










0.06 










0.01 


Jet corrections 






0.01 










0.01 


0.02 














Higher order terms 






0.01 










0.01 




0.01 






0.14 






Uncorrelated errors 




Bose-Einstein correlations 








0.01 










0.02 










0.01 




Colour reconnection 


















0.02 










0.01 




Background estimation 








0.01 










0.02 










0.01 




Monte Carlo statistics 






0.07 




0.24 






0.19 










0.17 




0.06 


Total 


0.01 


0.02 


0.20 


0.03 


0.25 


0.07 


0.16 


0.36 


0.08 


0.15 


0.01 


0.02 


0.27 


0.05 


0.07 



Table 7: Summary of systematic uncertainties for the couplings Agf , Ak 7 and A 7 . A description of the different sources is given in the text. Systematic 
uncertainties below 0.005 are indicated by a dash. 





Real 


Imaginary 


Source 






Hz 


\z 


7 
94 


9l 


9l 


9t 




Ay 


Hz 


\z 


7 
94 


9l 


9f 


9t 


Correlated errors 




Fragmentation 


0.04 


0.03 


0.01 




0.04 


0.01 


0.03 


0.05 


0.02 


0.02 


0.01 


0.01 


0.01 


0.07 




0.02 


W + W~ cross section 


0.02 


0.02 


0.01 


0.01 




0.02 






0.01 


0.01 


0.01 




0.04 


0.02 


0.02 


0.01 


Luminosity 


0.01 


0.01 






0.01 


0.02 














0.01 


0.02 


0.01 


0.01 


LEP energy 


0.01 








0.02 


0.04 


0.01 


0.03 


0.01 


0.01 






0.01 


0.02 




0.02 


W mass 












0.02 














0.01 


0.01 


0.01 




Calorimeter scale 


0.01 


0.01 


0.01 


0.01 




0.04 




0.03 


0.01 


0.01 






0.01 


0.01 




0.01 


Tracking 


0.02 


0.01 


0.01 


0.01 




0.02 






0.01 


0.01 






0.01 


0.02 


0.01 




Jet corrections 










0.01 


0.01 


0.01 


0.01 












0.01 




0.01 


Higher order terms 


































Uncorrelated errors 




Monte Carlo statistics 
Background estimation 


0.02 


0.02 


0.02 


0.01 


0.04 


0.04 


0.03 


0.02 


0.01 


0.01 


0.01 


0.01 


0.03 


0.05 


0.02 


0.03 


Total 


0.06 


0.05 


0.02 


0.02 


0.06 


0.08 


0.04 


0.07 


0.03 


0.02 


0.02 


0.01 


0.05 


0.09 


0.04 


0.04 



Table 8: Summary of systematic uncertainties for the combined evqq and \ivqq channels for C- or P-violating couplings. A description of the different 
sources is given in the text. Systematic uncertainties below 0.005 are indicated by a dash. 



Table 9: The combined results for 183 and 189 GeV for each W + W decay channel for the three couplings 
Ag^, Ak 7 and A 7 . The error includes the statistical and systematic uncertainty. 





Coupling 


Channel 






A 7 


evqq 


n nq+ - 09 


U - 4b -0.32 


n oi+O.ll 
°- 21 -o.io 


fivqq 


001 +o.io 

U ' Ui -0.10 


20+ - 64 
U ' 2U -0.34 


~ u - U8 -0.09 


rvqq 


M +0- 21 
U '^ i -0.37 


_ n ?1 +0.54 
u - a -0.39 


00 +0.18 


qqqq 


-0 03+ - 10 
U - Ud -0.10 


n O7+0.30 
u -^-0.26 


01+0.13 
u.ui_ 012 


Ivlv 


_n 17 +0.36 


_n 35+0-82 


U - U5 -0.13 



the systematic uncertainties. By this procedure the systematic uncertainties are folded with 
the proper statistical correlations between the optimal observables and the results of the 
fits include both the statistical and systematic errors. The systematic uncertainties listed in 
Table |7] are derived from the changes in the optimal observables mean values. They are not 
used as such in the analysis but serve only as a representation of the systematic contributions 
from the different sources. 

The systematic uncertainties for the maximum likelihood method are convoluted into the 
log L functions by assuming parabolic behaviour of the systematic errors around the fitted 
TGC value. 



6.3.3 Final results from W-pair production 

The combined results from all W + W~ decay channels at 183 and 189 GeV for the three 
couplings Agf , Ak 7 and A 7 , are obtained by combining the OO2 analysis of the evqq, \ivqq 
and qqqq final states with the OOi analysis of the rvqq and Ivtv final states. The correlation 
of the systematic errors between the different channels and energies are included as described 
in Section. |6.3.2| . The results for Agf , Ak 7 and A 7 , including systematic uncertainties, are 



listed in Table |9|. The final 68% and 95% combined W-pair result for the three couplings 
Ag^, Ak 7 and A 7 is summarised in Table |10[ The corresponding logL curves, including 
systematic uncertainties, are shown in Figure || 

The maximum likelihood method is used in the two semileptonic channels, evqq and 
fivqq, to measure the real and imaginary parts of the C- or P- violating couplings gj, g^, ky, 
and Ay, where V denotes either 7 or Z. The combined 183 and 189 GeV results, including 



the systematic uncertainties, for the C- or P-violating couplings are summarised in Table 10 



In all cases described above, each coupling is determined fixing the other couplings to 
their Standard Model values. The error intervals for each coupling are defined as the 68% 
confidence level intervals obtained by integration of the likelihood functions, to accommodate 
cases with non-parabolic behaviour of the log-likelihood function. 
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Table 10: Combined 183 and 189 GeV W-pair results for the three C- and P-conserving couplings, 
Ag^, Ak 7 and A 7 , and the C- or P- violating couplings. The error includes the statistical and systematic 
uncertainty. The corresponding 95% confidence level intervals are listed in the last column. 





fit result 


95% confidence limits 


As?? 


0.02+0-06 
—0.06 


[ -0.09, 0.14 ] 


A/t 7 


0.22+0-21 
—0.20 


[ -0.15, 0.66 ] 




0.04+0-06 
—0.06 


[ -0.08, 0.17 ] 


Rp(k ) 


-o lg+o- 1 ^ 

u - iy -0.17 


r n nisi 


Rp(\ ) 


17+0- 14 
u - i '-0.16 


[ -0 1 5 43 1 




-0 09+0-1? 
uuy — 0.11 


[ -0 SO 14 1 

1 U.OU ; U.14: 




07+0- 09 

u - u '-o.io 


r _n 19 n 95 1 


Rp(a"<) 


06+0-34 


\ -0 fi9 79 1 




-0 02+0-51 
uu ^— 0.51 


f -1 09 98 1 


Rp(n z ) 


n 07+9-?? 

u - u '-0.23 


1 u.oo, u .o\j \ 


Rp(n z ) 


—0 ofi+Q-?? 


U.U I ; U.JU 


111 tyrb'y J 


10+0-12 
u - iU -0.12 


r -0 14 o 33 1 


Tm(\ \ 

111 tyA'y J 


-0 08+0- 10 

u - U8 -o.io 


r n 97 n 1 1 l 


J 171 I K r? 1 

1 1 1 L\r\j^ J 


03+0- 08 


[-013 0191 


Im(\z) 


-0.03+0-07 
—0.07 


[ -0.16, 0.10 ] 


Im(gJ) 


n o 7 +0.30 
u - d '-0.30 


[-0.23, 0.95 ] 


Im(gl) 


-ooi+o- 57 


[ -1.10, 1.10 ] 


Im{gf) 


n 97 +0.20 
u -^-0.20 


[ -0.13, 0.65 ] 


Im(gi) 


0- 07 -oI 


[ -0.62, 0.76 ] 
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Table 11: Combined results for Ag^ , Ak 7 and A 7 from W + W~ production at 172-189 GeV, single-7 and 
single-W production at 161-189 GeV. The errors include systematic uncertainties. The corresponding 95% 
confidence level intervals are listed in the last column. 



Coupling 


fit result 


95% confidence limits 


Agf 


02 o+0.059 


[-0.087,0.141] 




n 022+0.119 
U.U22_q 115 


[-0.200,0.258] 


A*y 


U.U4U_q 052 


[-0.062,0.147] 



Table 12: Result of a three-parameter fit for Agf, Ak 7 and A 7 using the combined information from 
W-pair production at 172-189 GeV, single-7 and single-W production at 161-189 GeV. The statistical and 
systematic uncertainties are combined in a 68% one-dimensional error. The corresponding correlations are 
given in the last column. 



Coupling 


fit result 


Correlation 


Agf Ak 7 A 7 


Ag? 


O013+O.O66 
u.uid_oo 68 


1.0 -0.1 -0.6 
1.0 -0.1 
1.0 




04Q+0.110 
U.U43_ 11Q 


A-y 


023+ - 074 
U.U^_oq77 



7 Combined TGC results 

The measurements from single-7, single W and WW production, are combined with previous 
ALEPH results from W+W" production at 172 GeV 0, single-W production at 183 GeV [| 
and single-7 production at 183 [Q. The combined results are listed in Table |Tl|. In Figure || 
the corresponding one-parameter log L curves are shown. 

To study the full correlation between the parameters, two- and three-parameter fits, 
where two or all three couplings are allowed to vary, are also presented. The fits use the 
combined information from W-pair production, single-W production and single-7 production 
at 183-189 GeV. 

For the three parameter fit the results and the errors computed from a variation from the 
minimum of the logL functions of 0.5, are summarised in Table O including the systematic 



uncertainties. The correlation matrix of the three-parameter fit is also given in Table [12]. This 
correlation matrix is evaluated at the local minimum, and the correlations vary substantially 
depending on the exact value of the minimum. The projections onto the two dimensional 
plane of the three dimensional envelope of the 95% confidence level volume, representing the 
integration of the confidence over the corresponding third coupling, are shown in Figure |H| 
The 95% confidence limits of the respective 2-parameter fits of the three pairs of couplings 
(Ag^, A« 7 ), (Ag'f, A 7 ) and (Ak 7 , A 7 ) are shown as full lines. The systematic uncertainties 
are included in the limits shown. No deviations from the Standard Model expectations are 
observed. 
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8 Summary and conclusions 



The triple gauge-boson couplings have been measured using W-pair events at 183 and 
189 GeV, single-W production at 189 GeV and single-7 production at 189 GeV. Combining 
with previous ALEPH results from W + W~ production at 172 GeV, single-W production and 
single-7 production at 183 GeV, the three couplings Agf, Ak 7 and A 7 have been measured 
individually, assuming the two other couplings to be fixed at their Standard Model value. 
The results are 

-0.059 
-0.055 

0.119 
-0.115 

-0.054 
-0.052 ' 

where the error includes systematic uncertainties. The corresponding 95% confidence level 
limits, 



Agf = 


0.023 


Aft 7 = 


0.022 




0.040 



-0.087 < 


Agf 


< 


0.141, 


-0.200 < 


Ak 7 


< 


0.258, 


-0.062 < 


A 7 


< 


0.147, 



are in good agreement with the Standard Model expectation. Multi-parameter fits, where 
two or all three couplings are allowed to vary show also good agreement with the Standard 
Model. 

In addition, semileptonic W-pair events were used to set limits on the C- or P-violating 
couplings gj , ggf, Ry, and Ay, where V denotes either 7 or Z. No deviations from the Standard 
Model expectations are observed. 
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A Leptonic Neural Network Input Variables 



The Neural Network (NN) calculates an approximation of the multidimensional probability 
density function in the following 13 input variables, for signal and backgrounds. The NN 
is applied, after preselection, to events with at least two opposite charged tracks with 
momentum - after bremsstrahlung correction - in excess of 15 GeV. The NN uses variables 
related to the lepton candidates, to the missing momentum, global event variables, and WW 
kinematics. They are listed here together with their relative discriminating power, namely 
the statistical correlation with the neural network output: 

• missing mass squared (13.7%); 

• missing transverse momentum (11.5%); 

• angle between the two most energetic tracks (9.6%); 

• energy of the second most energetic track, (8.7%); 

• total energy found in a 12° cone around the beam axis (8.2%); 

• number of identified leptons with an energy greater than 15 GeV(8%); 

• missing transverse momentum with respect to the plane defined by the beam axis and 
the 3D-thrust axis (7.8%); 

• energy of the most energetic track (7.4%); 

• invariant mass of the two most energetic tracks (6.6%); 

• missing longitudinal momentum (6.3%); 

• scalar sum of the transverse components of the two most energetic tracks with respect 
to a 2D-thrust axis, built from the projection of the track momenta on the transverse 
plane (5.4%); 

• number of isolated neutral clusters with energy more than 4 GeV outside a cone of 10° 
around each of the two most energetic tracks and forming an invariant mass with each 
of them of more than 2 GeV (4.1%); 

• cosine of the angle between the most energetic track and the axis perpendicular to the 
plane defined by the second most energetic track and the z-axis (2.6%). 
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Figure 1: Distribution of a) the scaled energy, xe, and b) the absolute value of the cosine 
of the polar angle for single-7 events selected in 189 GeV data. The data are represented 
by solid dots, while the solid histogram shows the distribution for the Standard Model. The 
dashed and dotted histograms show the distribution for non-standard values of Ak 7 = ±5.0. 
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Figure 2: Negative log-likelihood curves, including systematic uncertainties, from the single- 
7 analysis for a) Ak 7 and b) A 7 for the 189 GeV data (dashed line), 161 - 183 GeV data || 
(dotted line), and the combined results (solid line). The curve for each coupling is obtained 
while fixing the other coupling to its Standard Model value. 
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Figure 3: The distribution of the lepton transverse momentum, p T , for single- W events 
passing the final selection cuts for the leptonic electron (upper left) and muon (upper right) 
W decay. The lower plot shows the visible mass distribution from single-W events passing 
the final selection cuts for the hadronic W decay. The data are represented by the closed 
circles. The histograms correspond to the Standard Model prediction. 
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Figure 4: The negative log-likelihood curves from the single-W analysis as functions of a) 
Ak 7 and b) A 7 for the 189 GeV data (dashed line), 161 - 183 GeV data [| (dotted line) 
and the combined results (solid line). The curve for each coupling is obtained while fixing 
the other coupling to its Standard Model value. Systematic errors are not included in these 
curves. 
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Figure 5: The distributions of the kinematic quantities cos 9 W , cos 9f, </>*, cos 9* et and 0* et from 
the combined sample of the evqq and fivqq channels at 189 GeV. The measured variables 
are the angle $w between the W~ and initial e~ in the W + W~ rest frame, the polar and 
azimuthal angles of the lepton, 9* and in the rest frame of its parent W, and the polar 
and azimuthal angles of a quark jet, 9? ct and 0*j e t, in the rest frame of its parent W. As no 
quark flavour tagging is performed each of the two ambiguous solutions enters with a weight 
of 0.5. The data are represented by solid dots, while the solid and dashed histograms show 
distributions for Standard Model and non-standard values of = ±0.5. 
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Figure 6: Distributions of the cosine of the W~ production angle, cosOw, at 189 GeV for 
a) ruqq, b) qqqq and c) iuiu events. The data are represented by solid dots, while the solid 
and dashed histograms show distributions for Standard Model and non-standard values of the 
TGCs. The shaded area represents the non-WW background. For qqqq events, each event 
enters with two solutions for cos 9w in the distribution with the weights P + and 1 — P + , 
where P + is the probability for a di-jet pair to be a W + . For iviv events, each event enters 
with two solutions for cos 6^ in the distribution with a weight of 0.5. 
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Figure 7: W-charge tagging distributions from qqqq W-pair events at 189 GeV. a) The 
distribution of the jet pair charge for W + (solid histogram) and W~ (dashed histogram) 
decays for Monte Carlo events, b) The probability P + as function of the charge difference 
between the two W's. c,d) Experimental distributions of Qw+ an d Qw- from semileptonic 
events. The data are represented by the dots and the Monte Carlo simulation by the 
histograms. The number of Monte Carlo events is normalised to the integrated luminosity 
of data. 
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Figure 8: The combined negative log-likelihood curves from the W-pair analysis of 183 and 
189 GeV data for the individual fits in the tvqq (dashed), qqqq (dotted) and tvtv (dashed- 
dotted) channels for the three couplings a) Agf, b) A/t 7 and c) A 7 . The curve for each 
coupling is obtained while fixing the other couplings to their Standard Model value. The 
systematic uncertainties are included. The combined result for all channels is shown as the 
solid curve. 
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Figure 9: The negative log-likelihood curves for the combined fits using single-7 (dotted), 
single- W (dashed-dotted) and W-pair (dashed) production at energies up to 189 GeV for 
the three couplings a) Agf , b) Ak 7 and c) A 7 . The curve for each coupling is obtained while 
fixing the other couplings to their Standard Model value. The systematic uncertainties are 
included. The combined result is shown as the solid curve. 
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Figure 10: Multi-parameter fits using the combined data from single-7, single-W and W-pair 
production at energies up to 189 GeV. The two-dimensional 95% confidence level contours 
for the three pairs of couplings, a) (Agf, Ak 7 ), b) (Agf, A 7 ) and c) (Ak 7 , A 7 ). The solid 
lines show the 95% confidence level contours of the two-parameter fit. The shaded area is 
a projection onto the two-dimensional plane of the three-dimensional envelope of the 95% 
confidence level volume. The Standard Model point is represented by a star. 
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